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Abstract 

The timing of the local group is used to test Modified Newtonian Dynamics 
(MOND). The result shows that the masses predicted by MOND are well below 
the baryonic contents of the Milky Way and the Andromeda galaxies. 
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1. Introduction 

Although dark matter has not been directly detected yet, there is strong evidence 
for its existence in the universe (Zwicky 1933, Rubin & Ford 1972). Recent 
reports showed that dark matter is essential in the formation of the large-scale 
structure of the universe (Springel, Frenk & White 2005). 

At the largest scale, the flatness of the observed universe and the cosmic 
microwave background (CMB) data requires 23% of the mass-energy density in 
the universe to be dark matter (Komatsu et al 2008). However, attempts to directly 
detect dark matter have failed despite decades of effort. 

It is reasonable to raise a question about dark matter being another kind of "ether" 
which was widely accepted at the end of 19 th century. The need for dark matter is 
based on Newton's law of gravitation and theory of general relativity, which are 
well-tested in the lab and up to the scale of the solar system. There is a possibility 
that both theories cannot be extended in some conditions. The most successful 
alternative to dark matter is Modified Newtonian Dynamics (MOND) as proposed 
by Milgrom(1983). 

MOND is extremely successful in explaining the discrepancy between the rotation 
and luminosity curves in spiral galaxies. In other words, MOND requires only the 
baryonic masses to explain the flat rotation curves of spiral galaxies. There is no 



need for the existence of dark matter in the galaxy scale (Sanders & McGaugh 
2002). Numerous observational data were used to test MOND in the last two 
decades (Pointecouteau & Silk 2005, Klypin & Prada 2009). Very few of them 
provide conclusive evidence to falsify MOND. The most convincing one is the 
recent gravitational lensing study of IE 0657-558 also known as the Bullet Cluster 
(Clowe et al 2006). The study shows that the location containing the most 
baryonic matter does not match the gravitational lensing centers, where most of the 
mass lies, after the collision of two galaxy clusters. This is viewed as convincing 
evidence for the existence of dark matter. 

In this paper I apply the timing of the local group to test MOND. I find that 
MOND predicts a much smaller mass for the Milky Way. The predicted mass is a 
factor of four smaller than the baryonic mass of the Milky Way. 

In Section 2, 1 repeat the calculation using the Newtonian motion and thus 
determine the timing between Andromeda (M3 1) and the Milky Way. I obtain the 
same results as in many previous papers (Kahn & Woltjer 1959; Mishra 1985; 
Binney & Tremaine 2008). 

In Section 3, MOND is used to calculate the motion and timing of these two 
galaxies. 

2. Newtonian 

The current separation between M3 1 and the Milky Way is about 740 Kpc (Binney 
& Tremaine 2008; Ribas et al 2005). They are moving toward each other with a 
speed of 125 Km s" 1 (Binney & Tremaine 2008). According to the Big Bang 
theory, these two galaxies started out close to each other, and then moved apart due 
to the expansion of the universe. The gravitational attraction slowed them down, 
halted, and then reversed their recessional motion. If we treat the Milky Way 
(with mass mi) and M3 1 (with mass m 2 ) as two point masses and assume they are 
moving in an almost radial orbit, their motion is determined by Newton's law of 
gravitation. If we choose a coordinate system where the origin is at the center of 
mass of the Milky Way and M3 1 system, and r x is the distance of the Milky Way 
to the origin, then 
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At r m (the maximum separation between M3 1 and the Milky Way), f=0. 

. 7 2GM 2GM ... 
r z = (2) 



r 4rdr 2GM r , 

By integration of the above equation we can find the time t 0m for M3 1 and the 
Milky Way to move from r = to r m . 
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Likewise, we can find the time t mn for M3 1 and the Milky Way to move from r m to 
the separation now r n ^=740 Kpcj. 

The total time for M3 1 and the Milky Way to move from to r = to r m and from r m 
to the current separation r n 



t mai = tom + t mn =13.7 Gyr = the age of the universe 
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From Equation 2, we can get r r 
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where r n = 740 Kpc is the separation between M3 1 and the Milky Way at the 
present time, and f n = 125 Km s" 1 is the corresponding speed. 

From Equations 3 & 4, we can find M numerically. 



M= m 1 + m 2 = 4.5X 10 1Z M 



This is far larger than the baryonic mass in these two galaxies which is in the order 
of 10 11 M (Binney & Tremaine 2008). This is considered as a strong evidence of 
the existence of dark matter. 



3. Modified Newtonian Dynamics (MOND) 

The original goal of MOND is to use only the observed matter to explain the flat 
rotation curves of spiral galaxies and the well-known Tully-Fisher relation. It is 
phenomenological rather than based on any sound physical hypothesis or principle. 
It formulates mathematical equations to fit the flat rotation curves of spiral galaxies. 
Regardless of this, MOND is quite successful and works well when applied to 
galaxies (Sanders 2008). 

If g m is the MOND gravitational acceleration, and g n is the Newtonian acceleration, 
Milgrom (Milgrom 1983; Sanders & McGaugh 2002) proposed that 

gm \l ( I gm I / O-0 ) = gn , (5) 

where ja(x)=l, when x> > 1 (the Newtonian limit), 
(x) = x, when x < < 1 (deep-MOND regime), 
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a = 1.2 X 10" cm s " is the MOND parameter which has the unit of acceleration 
(Sanders & McGaugh 2002). 

It was recognized that Equation 5 does not preserve the conservation of linear 
momentum (Felten 1984). In 1984 Bekenstein and Milgrom proposed a rigorous 
Lagrangian formulation of MOND. The modified Poisson equation becomes 

V • [ fi ( I V0 I /oq) V0 ] = 4ttG P , (6) 

where - V = g m . 

This modification conserves linear momentum and energy as it is derived from 
Lagrangian formulation. 



In deep-MOND regime, for high symmetry cases where the density distribution is 
spherical, cylindrical, or planar, Equation 6 reduces to Equation 5 which provides a 
much simpler algebraic relation between the MOND acceleration and the 
Newtonian acceleration. 

Thus, in the limit of low accelerations (deep-MOND regime), from Equation 5 the 
MOND acceleration is given by 

gm= V a 0§n 

The above simplified equation is a good approximation to compute the timing 
between M3 1 and the Milky Way, since most of the time the system is in the deep- 
MOND regime except for a very short period of time when these two galaxies 
were close to each other. 

The MOND equation of the motion between M3 1 and the Milky Way is given by 

d 2 r 
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where g n is the Newtonian acceleration given by Equation 1 , and r is the 
separation between M3 1 and the Milky Way. 
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This equation can be written as 
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At r m (the maximum separation between M31 and the Milky Way), f=0. 
f 2 = 2r\ y /a Gm 1 (lnr m — lnr) (7) 

f VlnJ-lnr = -^^^ f dt 

The left hand side = — r m yfn~ erf (^/lnr m — lnr ). 

The time for M3 1 and the Milky Way to move from r = to r m (the maximum 
separation) 



The time for M3 1 and the Milky Way to move from r m to the current separation r n 
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The total time for M3 1 and the Milky Way to move from r = to r m and from r m to 
the current separation r n 

ttotai = tom + 4« = 13.7Gyr 
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The maximum separation r m could be obtained by inserting the current separation 
r n = 740 Kpc and the speed f n = 125 Km s" 1 into Equation 7. 

f n 2 = 2n ja Gm 1 (lnr m -lnr n ) (9) 

From Equations 8 & 9, we can find the MOND mass m^f the Milky Way 
numerically. 

In order to find mjwe need to provide a value for 

^ /!+=*. 
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Assuming the mass of M31 (m 2 ) is about 1.5 times the mass of the Milky Way 
(mi) , we get 

m 1 = 1.2 X 10 10 M . 

This is only 24% of the estimated baryonic mass of the Milky Way which is about 
5X1O 1O M (Binney & Tremaine 2008) . 



Instead of using Equation 8 to calculate the masses, Equation 8 can be used to 
calculate the time required for M3 1 and the Milky Way to get to the current 
separation if the masses for the Milky Way and M3 1 are known. If only the 



baryonic masses of the Milky Way (5 X 10 1U M ) and M31 (1.5 X 5 X10 1U M ) 
are used, the time required for M3 1 and the Milky Way to reach to their current 
separation is 7.3 Gyr. Thus MOND's prediction of the age of the universe is much 
shorter than the believed age of the universe. This is understandable because the 
MOND force is stronger than the Newtonian force in the deep-MOND regime. It 
would take less time to complete the same motion. 

4. Discussion 

One of the uncertainties of the above argument is that we do not know the exact 
past history of the local group. For example, if one assumes that M3 1 and the 
Milky Way had already completed one orbit and are approaching the completion of 
the second orbit now, the MOND calculation shows that the total time is 17.5 Gyr. 
This is not consistent with the current age of the universe either. 

The other uncertainty is the assumption that M3 1 and the Milky Way are moving 
in an almost radial orbit. If the orbit is non-radial, the timing will be longer. In 
order to increase the timing from 7.3 Gyr (as predicted by MOND for the radial 
orbit) to 13.7 Gyrs, there should be a considerable amount of proper motion. The 
required proper motion will be in the same order of magnitude as the current radial 
velocity between M3 1 and the Galaxy. That is the proper motion in the order of 1 
Km s" . This is inconsistent with the current view of the universe that the initial 
motion between M3 1 and the Milky Way is due to Hubble expansion in the early 
universe which has zero angular momentum. There are only two massive 
members, M3 1 and the Milky Way, in the Local Group. It is hard to explain how 
the system acquired so much angular momentum. 
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It is interesting to point out that the proper motion of 10 Km s" is equivalent to 
10" 5 arcseconds position change per year of M31in the sky. This will be in the 
measurable range of SIM (Space Interferometry Mission) in the future. 



5. Conclusion 



MOND is extremely successful in explaining the flat rotation curves of spiral 
galaxies and the Tully-Fisher relation (Sanders & McGaugh 2002) without the 
need for dark matter. However, when it is applied to galaxy clusters, MOND 
predicts masses that are double (Sanders 2003) or even several times greater than 
the baryonic masses (Pointecouteau & Silk 2005). Sanders (2003) argued that 
MOND's prediction of more mass does not constitute the falsification of MOND 
because there might be other forms of matter in the system that are not visible 
(such as neutrinos). His argument is that if MOND predicts less mass than is 
observed, then it will be a definite falsification of MOND. More mass can always 
be found, but it is difficult to make observed mass disappear. 

In this paper I apply MOND to a simpler case, the motion and timing of M3 1 and 
the Milky Way. The result is that MOND predicts much less mass than the known 
baryonic mass for the Milky Way. The fact that MOND predicts less mass 
provides a potential problem for MOND. 

The manuscript originated from a term paper study of the "Galactic Dynamics" 
course at Rutgers University. I would like to thank Professor Sellwood for his 
invaluable comments as well as inspiration. I would also like to thank Rutgers 
University for letting me audit several astronomy courses in the past few years. I 
am grateful to Professor Stacy McGaugh for the interesting discussions on non- 
radial motion. 
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